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Interference of an array of independent Bose-Einstein condensates, whose experiment 
has been performed recently, is theoretically studied in detail. Even if the number of 
the atoms in each gas is kept finite and the phases of the gases are not well defined, 
interference fringes are observed on each snapshot. The statistics of the snapshot inter- 
ference patterns, i.e., the average fringe amplitudes and their fiuctuations (covariance) , 
are computed analytically, and concise formulas for their asymptotic values for long 
time of flight arc derived. Processes contributing to these quantities are clarified and 
the relationship with the description on the basis of the symmetry-breaking scenario is 
revealed. 

Keywords: Bosc-Einstein condensation; interference; Hanbury Brown and Twiss effect; 
snapshot; spontaneous symmetry breaking. 



1. Introduction 

Interference fringes are observed when two independently prepared Bose-Einstein 
condensates (BECs) are released and overlap after a time of flight.^ The simplest 
explanation of this phenomenon is based on the spontaneous symmetry breaking 
of BEC, where it is assumed that the phases of the two Bose gases are individually 
fixed upon condensation, which enable the gases to interfere with each other.12131 

When the number of the bosonic atoms in each gas is finite, however, the sym- 
metry of the system can never be broken. Javanainen and Yoo pointed out in their 
seminal pape?^ that, even if the number of the atoms in each gas is well defined 
and its phase is accordingly uncertain, the two gases can exhibit interference. Note 
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that the quantum mechanics predicts rcsuh of accumulation of many experiments. 
We usuaUy look at the single-particle distribution to study interference in the den- 
sity profile, which predicts the image obtained by accumulating many photos of the 
overlapping BECs. No interference can be observed in this quantity if the number of 
the atoms in each gas is fixed and its phase is not well defined. However, Javanainen 
and Yoo demonstrated that interference fringes show up on each snapshotW It is 
important to notice that we look at many identical atoms at the same time by 
taking a photo of the cloud. The symmetrization of the wave function gives rise 
to correlations among the multiple atoms in the cloud (Hanbury Brown and Twiss 
effectP^), which make the probability to find an atom in the presence of other iden- 
tical atoms non-unform in space. This is the origin of the oscillation in the den- 
sity profile o il each sn apshot . Such an interference is called "measurement-induced 
interference"l2EinEIHI gince measurement (taking a photo) extracts an interference 
pattern. 

An experiment was carried out to study the interference of an array of indepen- 
dent BECs.^Note that the snapshot interference patterns differ from shot to shot0 
A certain number of snapshots were collected by repeating the experiment, and the 
fluctuations of the fringe amplitudes were analyzed, on the basis of the hypothesis 
of the broken symmetry of the BECs. The interference of the array of BECs is 
briefly discussed in Ref. 1101 concentrating on the fluctuation of a measure of inter- 
ference, with each BEC consisting of a fixed number of atoms, without assuming 
the symmetry breaking. In this paper, we set up tools for studying the snapshot 
interference. These are essentially equivalent to the ones suggested in Refs. [TOland 
[TTl where the interference of the array of BECs, each containing a fixed number of 
atoms, is studied. The time-evolution of the fringe amplitudes expected on snap- 
shots and their fluctuations (covariance) are computed, and concise formulas for 
the average fringe amplitudes and their covariance are derived, which characterize 
the statistics of the snapshot interference patterns after long time of flight. The 
processes contributing to these quantities are clarified and the relationship with 
the description on the basis of the symmetry-breaking scenario, which well explains 
the experimental data, is revealed. 



2. Snapshot Profile and Its Fluctuation 

Suppose that there are N identical bosonic atoms and one takes a "photo" of 
the cloud. The probability of finding the N atoms at positions {ri, . . . ,rN} at an 
instant t is given by the A^-particle probability distribution function 

Pi^^(ri, . . . , rjv) = -^ii^Hri) ■ ■ ■ V'^(r^)^(r^) • • • V-Cn))*, (1) 

''In particular, the fringe pattern shifts from run to run, and the fringes are smeared out if all the 
snapshots are superimposed, reproducing the quantum-mechanical expectation for single-particle 
distribution. 
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where i/^if) is the field operator of the bosonic atom, satisfymg the canonical com- 
mutation relations 

mr),^:\r')]^S\r^r'), etc., (2) 

and the average (■••)* is taken over the state of the cloud evolved from time to 
t, while the operators are time-independent (Schrodingcr picture). This probability 
is normalized to unity, 

" J3^. ... j3^., pW 



d^ri...dVArPr'(n,...,rjv) = 1. (3) 



Note the relation 



J d'reP,^''\n,...,n,...rN) = P^^''-'\n,...,n-i,n+,,...rN). (4) 

Given a single configuration {ri, . . . , rjv}, the snapshot density profile of the cloud 
is constructed as 

1 ^ 



N 

i=l 



This profile is also normalized as 



J d^rp{r) = 1. (6) 

Notice that the specific configuration {ri,...,rjv} appears according to the 
probability P^ ^\ri, . . . ,rN) given in ([T]) and the snapshot profile p{r) fluctuates 
from run to run. The average profile 

W)^ Jd^r,---d^rMPt^''\r,,...,rN)pir) = P,^'\r), (7) 

which represents the image obtained by accumulating many snapshots, coincides 
with the single-particle probability distribution pj^\r) [We omit to specify the 
time dependence of the averages • • • , which depend on the time specified in Pj^^'' ; 
p{r) itself is time-independent]. However, p{r) is not the quantity of interest. We 
arc interested in the presence of interference fringes on each snapshot p{r)^ while 
the single-particle distribution P^^\r) (superposition of many snapshots) does not 
exhibit interference pattern.SI In order to detect an oscillating pattern on each 
snapshot, we look at the Fourier spectrum of the density profile p{r), 

p{k)^ j d^rp{r)e-'^ ''. (8) 

Possible nontrivial spikes in p{k)^ representing the presence of interference fringes 
in each snapshot would disappear if average is taken over all the snapshots. The 
reason for the disappearance of the interference fringes is the random shift of the 
interference pattern from snapshot to snapshot, which smears out the fringes .SI We 
can discard such random offset by removing the phase of the Fourier spectrum, i.e.. 
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by looking at the (square) modulus The fringe spikes would then survive 

on average in the quantity 



Stik) = |p(fc)|2, (9) 

if the majority of the snapshots exhibit interference patterns with a definite fringe 
spacing, irrespectively of the random spatial shifts. The fluctuation of the fringe 
spectrum from snapshot to snapshot is estimated by the variance, or more generally, 
by the covariance 



Ct{k, k') = |p(fc)|2|p(fe')|2 - \~p{k)? ■ \~P{kT- (10) 

These are our tools for studying the snapshot interference patterns. Notice that we 
could define an observable representing a (relative) phase 9{k) = argp(fe), which is 
not conjugated to the (relative) density. 
By noting 

1 ^ 

-^w-^^E^^'""' (11) 

1=1 

one realizes that these quantities are controlled by few-particle distribution func- 
tions. Indeed, one gets 

Sm = ^I?\k) + j^. (12) 
- W(N^fi''^^^^'^ ~ ^-^^l[^\k)l[^\k') + 0{j,) , (13) 

where 

l\^\k)= y dVid^r-jPf )(ri,r2)e*'=(''i-''=), (14) 

4^\k,k') = J £riSr2d^r:id^riPi^\ri,r2,r3,ri)e'^-^'-^-'-^-^+''''<''-'-'-^\ (15) 

The interference pattern of the N particles is essentially ruled by the two-particle 
distribution P^'^\rl^r2)^ while its fiuctuation by p/"'''(ri, r2, r3, r4). The formulas 
and show that the present formulation is essentially equivalent to the one 
employed in Ref. [TTl when the number N of the atoms in the cloud is large. 

Note the symmetries of these quantities, St{k) = St{—k), Ct{k,k') = 
Ct{k',k) ^ Ct{k,-k') ^ Ct{-k,k') ^ Ct{-k,-k'). 



3. Interference of an Array of Bose-Einstein Condensates 

Let us now study the interference of an array of independent BECs. Suppose that K 
independent BECs are created at regular distances in a periodic lattice with a lattice 
constant d, forming a ID array of BECs along the x axis (in the actual experiment, 
the array is loosely confined in the transverse directions and the BECs are like 
disks). We assume that each BEC contains exactly N atoms of mass m, which are 
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all condensed in the ground state of each well at zero temperature. We consider 
ideal gases, neglecting the intra-atomic interactions. The wells are well separated 
with vanishing tunneling probability between adjacent wells. We emphasize that 
these BECs are mutually independent, with the fixed number of atoms in each 
BEC endowed with no definite (relative) phase, and there is no phase correlation 
among them. At time t = 0, the optical lattice is switched off and the array of 
BECs is released. After free expansion in space (we do not consider gravitational 
field), a photo of the cloud is taken to observe an interference pattern among the 
BECs. The expansions of each BEC in the longitudinal and transverse directions 
are decoupled, and in addition, the expansion in the transverse directions is slower 
than that in the longitudinal direction, since the transverse confinement is weaker 
than in the longitudinal direction. In the following analysis, we concentrate on the 
dynamics in the longitudinal direction, which is relevant to the appearance of the 
interference fringes. 

Let ipn{x) {n — 1, • . ■ , K) denote the wave function of the ground state of the 
nth potential well and a„ the associated annihilation operator. We assume that the 
wave functions are orthogonal to each other. 



dx (^* {x)ipn' (x) = Snn' , (16) 

and the bosonic operators of different wells commute, 

[a„,a,\,] = 6nn', etc. (17) 
During the time of flight, the wave functions evolve in free space according to 

(Pnix,t) = e^'i^(p„{x). (18) 



Let us first look at the average of the snapshot profiles, p{x), which gives the 
single-particle distribution of the atoms in the cloud, ^'(x), as shown in ([7]). It is 
given by 

n=l 

in which no interference fringes are observed. This is the definition of the "indepen- 
dence" of the BECs. Interference patterns are however found on each snapshot. Let 
us look at the average spectrum St{k) and the fluctuation C((fc, k') of the snapshot 
profiles, introduced in (|9))- p0)) and evaluated by (fT2|) - ((T3| . Since the two-particle 
distribution of the atoms after the time of flight is given by 



KN 



KN - 



I 1 ^ 

- P[^\x{)P['-\x2) - ^ E l¥'n(xi,i)|Vn(^2,i)l' 
\ n=\ 

^ K K \ 

+ E E V',*2(2^i,0v^i(a^2,<)<y5n2(a^2,<)V'ni(xi,<) 1 , 

"1 = 1 "2 = 1 / 

(20) 
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its Fourier transformation (jl4p yields 



1 ^ 1 ^ ^ \ 

n—1 ni — 1 712 — 1 / 



n—l n 

ni/n2 

(21) 

where 

1 ^ r 

n=l 

are the Fourier transforms of the single-particle distribution p}^\t) and of the so- 
called "interference term," respectively. lj:'^\k,k') defined in is also composed 
of Xijik, t) as 

^* ^'"'^^^ (i^TV)! 

X a/ 7V(iV - (5„i„J(A^ - (S„i„3 - (5„2„3)(A^ - i5„i„4 - (5„2n4 - <5„3„J 

{k',t). (23) 

3.1. Gaussian Wave Packets 

Let us look explicitly at the time-evolution of the interference of the array of BECs. 
In the experiment reported in Ref. [HI -fsT = 30 condensates of ^^Rb, each containing 
N 10* atoms, were created with a lattice constant d ~ 2.7 /im. Each condensate 
can be regarded as being confined in a harmonic potential with frequency a;/27r = 
4 kHz in the longitudinal direction (while with ujj_/2'k = 74 Hz in the transverse 
directions), and its wave function in the longitudinal direction x would be well 
approximated by a Gaussian 

^f^) = 1 ,-(x-nrf)V4(A.)^ (n = 1, . . . , if ) (24) 

' ^27r(A.T)2 ^ ' ' ^ \ J 



of size Ax = \jyiJ^2Amjj = 120 nm. The overlap between the Gaussian wave functions 
of diff'erent wells, 

" dx (x)^„, (x) = e-("-"^)''^V8(A.)^ ^ ^25) 

is negligibly small since Ax ^ d, and they are approximately orthogonal to each 
other [see (HH)]. 

For these Gaussian wave functions, the Fourier transform of the interference 
term Xn,ini{k^t) defined in is given by 

Xn n {k t) = g-('^^)^fcV2g-"^t^(Ax)2[fc-(ni-n2)md/rit]V2g-i(ni-|-n2)M/2 ^2Q) 
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and by inserting it to l[^\k) in (HI]), wc get the average spectrum St{k) at an 
instant t, 



St{k) 



\K^siT?{kd/2) KN 



K-l 



where 



md 



1 



ht 



(27) 
(28) 



This expression is exact and vahd for any and t. 

The time-evolution of the average spectrum St{k) in ([27)) is shown in Fig. [T] 
with the experimental parameters recapitulated above. After a short transient 
period, the average spectrum St{k) exhibits 2K — 1 = 59 peaks. The central peak 
at fc = represents the overall profile of the whole cloud, which becomes narrower 
and narrower as the cloud expands during the time of fiight, while its height is kept 
constant ^((O) = 1, since the density profile p{x) is normalized to unity in ([6|)f^ 
The side peaks at fc ~ ±fc„ {n = 1, . . . , K ~ 1), on the other hand, represent the 
interference fringes in the snapshot profiles. Those peaks shift toward the center. 



^St{k) given in II27II for the Gaussian wave packets yields St{0) ^ 1, contradicting with the 
normaUzation condition S't(O) = 1 in This is due to the nonorthogonality of the Gaussian 
wave functions ipn{x) in l|25| l. The deviation of ^((O) from unity is actually of order of the overlaps 
among the Gaussian wave functions and can be neglected in the regime Ax <C d assumed in the 
present analysis. 




l.O'O 



Fig. 1. The time-evolution of the average spectrum St{k) of the interference of an array of 
Gaussian BECs, with the experimental parameters K = 30, A'^ = lO**, d = 2.7 ^m, Ax = 120 nm, 
uj/2n = 4 kHz. 
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Fig. 2. The covariance Ct{k, k') of the spectrum of the interference of an array of Gaussian BECs 
at time t = 22 ms, with the same experimental parameters as in Fig. [T] 

according to the formula for fc„ in (j28p . while they become sharper and higher as 
time goes on. In the experiment reported in Ref. [9l snapshots of the density profiles 
were taken at t = 22 ms after the release from the optical lattice. The covariance 
Ct{k, fc') of the interference spectrum at this time is shown in Fig. El for a large TV. 



3.2. Asymptotic Behavior 

The asymptotic heights of the peaks in St{k) (Fig. [1]) and those in Ct(fc, fc) (Fig. [2]) 
for large t can be estimated without assuming the Gaussian wave packets. Observe 
the asymptotic behavior of the wave packet for large t, 

<fn{x,t)= J dfce-*'''='*/2mgife(a;-nd)^(^) 

_ ^^^^mi.^ndr-/2M~ (11^^.^ nd)) . (29) 

V iht \ht I 

Then, the Fourier transform of the interference term defined in (|22l) . 

Xnxn.'i^it) ^ g-i{ril-nl)md^ /2ht 

-i{ht /m)[k—{ni— 712 )rnd/ht] k' 



X ldk'r{k'-n,'^]Jk'-n,^ 



e 

(30) 



becomes sharply peaked at fc = (ni — n2)md/ht with unit height Xnin2(^i^) ~ 
^-i(nl-nl)md' /2M ^ and l[^\k) givcn in ^ develops 2K - 1 peaks 

' 1 at fc = 0, 

ll'\k)-{ K-n , , (31) 



I 



I 
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where fc„ is defined in (|28|1 . At these wave numbers, l['^\k,k') given in is 
estimated to be 



iV4(ifiV-4)! 



{KN)l 



N'^iKN - 4)! 



(K ~n)[ 2K - 2n - 1 



(m = 0,n > 0), 
6 1 \ 4n" 

(m = n > 0), 

(A' - m){K -n) + 2{K ~m~ n) (^1 - 1^ - ^ 

(m > n > 0), 
(32) 

for large t. Therefore, combining these expressions, we get concise expressions for 
the peaks in the average spectrum St{k) and their fluctuations Ct(fc, k') for large t, 



(A'iV)! 



Stikn) 



K 1 



KN 



(n = 0), 
{n > 0), 



(33) 



Ctikmi kn) ^ < 







iK~n){K-n-l) f 1 



2(K-m-n) ^( \ 



(?7i = 0,n > 0), 
(to ~ n > 0), 

{m> n> 0). 



(34) 



Remember the symmetry Ct{k, k') = Ct{k', k). 

The peaks of St{k) at fc 7^ stem from the "interference terms" Xnin^i^jt) 
in the last term of l[^\k) in ((2T|) . This represents the interference of the atoms 
originating from the nith and n2th condensates, and it exhibits a peak at fc = 
{ni—n2)md/ht, as shown in (pO)) [see also (|26p ]. This is because the atoms interfering 
between these condensates at time t are those which have travelled the distance jni — 
n2\d/2 in time t, i.e., they have velocities =f("-i —n2)d/2t, and accordingly momenta 
=p(ni — n2)md/2t. The interference between the atoms counter-propagating with 
these momenta exhibits oscillation with wave number k = (ni — n2)md/ht. There 
are K — 1 different distances between condensates, and hence there are 2{K — 1) 
nontrivial peaks representing the interference fringefEl plus the peak at fc = 0. 
Therefore, there are 2K — 1 peaks in the spectrum S't(fc). 

The peak at k = kn represents the interference of particles originating from pairs 
of condensates which are separated by a distance nd. There are K — n such pairs of 
sources contributing this interference, among which there is no phase correlation. 



•^Notc that the density profile p{x) is a real function and its Fourier spectrum |p(fc)p is an even 
function of fc. 



18, 2013 22:55 WSPC/INSTRUCTION FILE BECarray 



10 S. Ando, K. Yuasa and M. lazzi 



^» St{ki) 



■• 

■I 



15 

K 



Fig. 3. The asymptotic strength of the average spectrum St(k) at fc = fci and its fluctuation 
Ct{ki, fci), as functions of the number of condensates K, in the large N limit. 



That is why the spectrum S't(fc„) is proportional to K — n, without interference 
among the contributions from different pairs of condensates, while its denominator 
is due to the normahzation of the density profile. On the other hand, the 
fluctuation Ct{k,k') is controlled by the number of independent ways of selecting 
two pairs of condensates separated by a distance nd, which is given by (K — n) [K — 
n-l). 

In Fig. 131 the asymptotic strength of the average spectrum St{k) at fc = fci 
(interference between adjacent condensates) and its fluctuation Ct{ki, ki) are shown 
as functions of the number of condensates K. Both average spectrum St{ki) and 
fluctuation Ct(fci, ki) decay monotonically as K increases, except for the fluctuation 
Ct{ki, fci) = for K ~ 2. For the double-condensate case K = 2, the interference 
pattern with perfect fringe contrast, which corresponds to St{k) = 1/ 4, is certainly 
observed on every snapshot with vanishing fluctuation Ct(fci, fci) = 0.^^ 



3.3. Spontaneous Symmetry Breaking 

The formulas for (/fc) and 4^\k,k') in ^ and and hence St{k) and 

Ct(k, k'), can be reproduced on the basis of the "symmetry-breaking scenario" em- 
ployed in Ref. [HI to analyze the experiment. This explains well the data. Suppose 
that the U(l) symmetry of the system is spontaneously broken upon condensa- 
tion and each condensate is described by a "macroscopic wave function" ipn{x)e^^"^ 
{n ~ 1, • . • , K), where 0„ is the phase of the condensate, which is randomly fixed 
upon condensation, varying randomly from to 27r from one realization to another. 
Note that in this scenario the number of atoms in each BEG is not fixed. The snap- 
shot profile for a single realization of the array of condensates, with a given set of 
{01, ... , 9k}, is then written as 

2 

(35) 



K 
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and the average spectrum of many snapshots for many different reahzations of the 
array of condensates is computed as 



K K 

= \It'\k)\^ + ]^ E E (36) 

ni — 1 n2— 1 

where lj:^\k) and Xnin^ik^t) are defined in p2|). This coincides with ([21]) and (|27p 
in the large N hmit. It is also the case for /j^' (fc, k') and Ct(fc, fc'). The average over 
the phases induces the constraints on rii and n[ as in (|23p. and the formula psp is 
recovered once 1 and iV + 2 are approximated by N in the large N limit. In this 
way, the description on the basis of the spontaneous symmetry breaking provides 
a good approximation in the large TV limit.^ 

From this point of view, the decay of the strength of the average spectrum 
St{kn) in ([33]) and that of the fluctuation Ct{km, kn) in (|34|) as functions of K are 
understood in the following way. If the number of the condensates K is large, it 
is very rare that the condensates are endowed with an appropriate set of phases 
{^1, . . . ,0k} which allow a constructive interference among the condensates. That 
is why the average spectrum St{kn) becomes smaller for larger K, as shown in Fig. 
[3l At the same time, snapshot patterns with less contrast become more typical 
for larger K , and the fluctuations of the fringe spectrum Ct{km, kn) is accordingly 
reduced as K is increased [while its ratio to the average spectrum S't(fc„) remains 



4. Summary 

We have discussed the interference of an array of independent BECs, each consisting 
of a fixed number of atoms, without assuming the symmetry breaking. The statistics 
of the snapshot interference patterns, i.e., the average fringe spectrum expected on 
each snapshot and its fluctuation (covariance), have been investigated. An analytical 
formula for the average spectrum valid for any time and any number of atoms has 
been presented]^ and concise formulas for both average spectrum and covariance 
valid for large time have been derived. The processes contributing to these quantities 
have been clarified, and the formulas have been shown to coincide with the ones 
based on the symmetry- breaking scenario if the number of the atoms is large enough, 
which explains well the experimental data. 

An interesting aspect of the present subject is to discuss snapshots, while quan- 
tum mechanics predicts result of accumulation of many experiments. As pointed 
out in Ref. llOi the average spectrum analyzed in the present paper does not supply 
definite information about each single snapshot, since the variance is not vanishing. 



"^An analytical formula for the covariance is also available but is not presented here for brevity. 
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It would be a challenging future subject to explore a better strategy to study the 
snapshot interference in a more conclusive way. 

Another interesting issue would be how it is possible to experimentally discrim- 
inate the two scenarios, i.e., measurement-induced interference and spontaneously 
symmetry breaking, and how they are related to each other. In principle, the dif- 
ference would become prominent when the number of atoms N is small, but at the 
same time, the visibility of interference pattern would be low with small N and 
experimental observation would be difficult. It would also be a challenging future 
subject to seek for a strategy to clarify the difference between the two scenarios. 
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